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ABSTRACT  OF  REPORT 


Backscattering  from  Three 
Dimensional  Conducting  Bodies 

by 

Arunkumar  P.  Nayak 
and 

Hagop  Panossian 

University  of  California,  Los  Angeles 
School  of  Engineering  and  Applied  Science 

Principal  Investigator 
Professor  Cornelius  T.  Leondes 


It  is  our  endeavor  here  to  develop  a  more  reliable  and 
complete  Synthetic  Aperture  Radar  target  modeling  capability.  The 
algorithm  to  solve  such  problems  should  be  practically  feasible 
and  accurate.  It  must  also  have  wide  range  of  application. 

The  main  core  of  the  study  is  to  calculate  the  electro¬ 
magnetic  (EM)  backscattering  from  life-size  objects  like  a  car  or 
a  tank,  etc.  To  accomplish  that,  we  will  pursue  high  frequency 
methods. 

The  three  most  popular  high  frequency  methods  are  discussed. 
They  are  (i)  Physical  Optics  together  with  Physical  Theory  of 


Diffraction,  (ii)  Numerical  methods,  and  (iii)  Geometrical  Optics 
in  conjunction  with  Geometrical  Theory  of  Diffraction.  We  found 
that  Geometrical  Optics  together  with  Geometrical  Theory  of 
Diffraction  is  the  most  practical  way  of  calculating  EM  Back- 
scatter  from  conducting  objects.  Also,  it  is  easier  to  apply  and 
has  conceptual  simplicity. 

After  making  the  choice  of  the  method,  we  applied  them 
to  calculate  EM  backscatter  from  a  sphere  and  a  spherical  segment. 
The  latter  object  is  a  simplified  model  of  a  car. 

We  have  given  all  the  equations  in  detail  for  three  cases 
and  also  we  have  given  Hierarchical  plus  Input-Process-Output  charts 
for  the  case  of  overview  and  coding  purposes. 


CHAPTER  ONE 


INTRODUCTION 

It  is  our  endeavor  here  to  develop  a  more  reliable  and  complete 
Synthetic  Aperture  Radar  target  modeling  capability.  The  algorithm  to 
solve  such  problems  should  be  practically  feasible  and  accurate.  It 
must  also  have  wide  range  of  application. 

The  main  objective  of  the  study  is  to  calculate  the  electro- 
magnatic  (EM)  scattering  from  life  size  objects  like  a  car  or  a  tank, 
etc.  To  accomplish  that  we  will  pursue  high  frequency  techniques. 

When  the  wavelength  is  small  compared  to  the  characteristic  dimensions 
of  the  scattering  body,  which  is  so  in  our  case,  asymptotic  high  fre¬ 
quency  methods  must  be  used. 

We  begin  with  the  basics  of  the  electromagnetics,  namely  the 
Maxwell's  equations  and  move  on  to  the  discussion  of  the  three  most 
popular  high  frequency  methods.  These  three  methods  are  (i)  Physical 
Optics  together  with  Physical  Theory  of  Diffraction,  (ii)  Numerical 
Methods,  (iii)  Geometrical  Optics  in  conjunction  with  Geometrical  Theory 
of  Diffraction.  After  doing  some  deep  probing  we  have  come  to  a  con¬ 
clusion  that  Geometrical  Optics  combined  with  Geometrical  Theory  of 
Diffraction  is  the  most  practical  method  for  calculating  EM  scattering 
from  objects.  The  method  is  general  enough  to  consider  any  kind  of 
polarization. 

After  making  a  choice  of  the  method  we  applied  it  to  calculate 
the  EM  backscattering  from  a  sphere  and  a  spherical  segment.  The 
latter  object  is  a  simplified  model  of  a  car. 


CHAPTER  TWO 


THEORY  OF  SCATTERING 

The  electromagnetic  field  at  a  time  t  and  at  any  point  in  a 
linear,  homogeneous  and  isotropic  medium  of  electric  permittivity  e, 
magnetic  permeability  u,  and  zero  conductivity  is  described  by  the 
homogeneous  Maxwell's  equations 


(1) 

* 

7  *  E  =- y  cFt 

(2) 

and  auxiliary  equation 

7  *  H  =  7  •  t  =  0  (3) 

which  govern  the  behavior  of  the  electric  field  ^  and  of  the  magnetic 
field  H  at  all  ordinary  points  in  the  space. 

In  the  following  we  shall  consider  only  the  particular  case  of 
monochromatic  radiation. 

If  a  scattering  body  is  embedded  in  the  medium,  (1),  (2),  and 
(3)  are  satisfied  by  the  incident  or  primary  fields  l:1  and  H1,  by  the 
total  (incident  plus  scattered)  fields  E  and  H  and  therefore  also  by 
the  secondary  or  scattered  fields  1:S  and  Hs,  which  represent  the  dis¬ 
turbance  added  to  the  primary  fields  by  the  scatterer. 

The  wave  number  k  and  the  intrinsic  impedence  Z  of  the  medium 
surrounding  the  scatterer  are  given  by 


Z  =  Y*1  =  /H7c 


(5) 


i 

& i 


l» 

J 


where  w  is  the  angular  frequency  ana  X  the  wavelength.  The  time 
dependence  factor  exp(-jut)  is  suppressed  throughout.  So  the  Eq.  (1) 
and  (2)  will  become 

7  x  h  =  -jkYi:  (6) 

7  x  £  =  jkZH  (7) 


2.1  Wave  Propagation 

For  harmonic  time  dependence  it  can  be  shown  from  Eq.  (1), 
(2),  and  (3)  that 


*  *  -*• 

(V  ♦  k  )H  «  0 


(7**  k*)£  =  0 


with  the  auxiliary  condition 


7  •  H  =  7  •  £  =  0 


2.2  Boundary  and  Radiation  Condition 

We  are  going  to  limit  our  considerations  to  perfect  conductors, 
i.e.,  we  require  that  the  tangential  component  of  the  total  electri- 
field  at  any  regular  point  of  the  scattering  surface  be  zero. 

t  -  (£  •  n)n  =  0  (10) 


where  n  is  the  unit  normal  to  the  surface  directed  from  the  body  into 
the  surrounding  medium.  It  is  seen  from  Maxwell's  equations  that 


condition  (10)  implies  that  the  normal  component  of  the  total  magnetic 
field  is  zero  at  any  regular  point  of  the  surface  scatterer. 


2.3  Integral  Equations 

For  an  electromagnetic  wave  incident  on  a  body,  integral  equa¬ 
tions  can  be  derived  from  which  to  determine  the  fields  induced  on  the 
surface  of  the  scatterer. 

It  is  convenient  to  confine  attention  to  a  three-dimensional, 
closed  and  bounded  body  whose  surface  S  is  regular. 

Stratton  [1]  has  given  the  necessary  equation 

2(r)  =  l:1  (r)  +  J  {jhZCn^fi)^  +  (n^)  x  Vip  +  ds  (11) 

S 

H(r)  *  HV)  -  Jp  J{ jkYfrjXE)*  -  (n^H)  x  V*  -  (^ -il)V^}ds 

S 

(12) 

Maue  [2]  has  given  an  expression  for  surface  current  which  follows 
from  Eq.  (12) 

2n(r)  xfiV)  =  3(r)  -  n(r)x  jftr^  x  dSj  (12a) 

S 


4>  =  exp(jkR)/R  (13) 

R  =  |r  -  (14) 

r  r.  e  S 

f  * 


The  del  operator  operates  on  the  surface  pt.  r^  only. 


If  the  body  is  perfectly  conducting,  the  surface  current 

density 

J  =  2n  x  H1  (15) 

and  this  follows  immediately  from  Eq.  (12a). 

2.4  Physical  Optics  (PO) 

"Physical  Optics"  denotes  an  approximate  method  for  the  deter¬ 
mination  of  the  field  scattered  by  an  object  through  an  assumption 
about  the  specific  form  of  the  field  distribution  on  the  surface. 
Explicitly,  it  is  a3sumed  that  the  field  there  is  the  geometrical 
optics  surface  field,  implying  that  at  each  point  on  the  geometrically 
illuminated  side  of  the  body  the  scattering  at  the  surface  takes  place 
as  though  from  an  infinite  tangent  plane  at  the  point,  whereas  over 
the  shadowed  portion  of  the  body  the  surface  field  is  zero.  For  a 
perfectly  conducting  body  the  postulated  current  distribution  is 
therefore 

J  *  2n  x  ft1  on  the  illuminated  side  (15) 

«  0  on  the  shadowed  side. 

A 

where  n  is  a  unit  vector  normal  drawn  outwards  as  regards  body.  H 
is  the  incident  magnetic  field.  The  scattered  fields  radiated  by  these 
currents  are  given  by  Rusch  and  Potter  [3], 

£(p)  *  -  ^  \  tf*V)V<J>}ds  (16) 

S  k 
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(17) 


ft(p)  =  4T  X  V^ds 

S 

when  r  »  X  we  can  approximate 

’ 

ES(p)  =  -  4>  [J  -  (J*aR)aR]  exp  (jkrj*aR)ds  (18) 

S 

fiS(p)  =  ii-  ♦  [^*aRl  exP  (jkr*aR)ds  (19) 

< 

s 

jp  =  exp  (- jkR)/R  ,  R  =  |r[  (20) 

In  contrast  to  the  geometrical  optics  expression  for  the 
scattered  field,  the  above  integrals  are  freauency  dependent,  and  it 
is  therefore  possible  that  physical  optics  provides  a  more  accurate 
estimate  of  the  scattering. 

Physical  optics  provide  an  accurate  estimate  of  scattering 
whenever  this  arises  from  a  portion  of  the  surface  where  the  actual 
surface  field  is  adequately  approximated  by  the  postulated  distribution. 
If  a  body  has  one  (or  both)  of  its  radii  of  curvature  infinite  at  the 
specular  point,  the  geometrical  optics  cross  section  is  infinite,  and  a 
particular  advantage  of  the  physical  optics  is  that  a  bounded  (and 
wavelength  dependent)  estimate  of  the  scattering  cross  section  can  be 
now  estimated. 

2.5  Physical  Theory  of  Pi  'fraction  (PTD) 

PTD  is  an  integrative  technique  in  which  the  surface  field  of 
Geometrical  Optics  (GO)  is  refined  by  the  addition  of  a  nonuniform  edge 
component.  The  method  is  therefore  an  extension  of  Physical  Optics  (PO) 


and  since  the  nonuniform  component  is,  in  principle,  approximated  in 
terms  of  the  surface  field  on  a  (local)  wedge,  the  result  should  con¬ 
stitute  a  caustically  corrected  version  of  the  GTD  prediction.  A 
description  of  the  method  is  given  in  Ufimtsev* s  monograph  [4]. 

Ufimtsev  postulates  that  the  total  induced  surface  current  J 
can  be  represented  as  the  sum  of  two  components,  one  the  physical 
optics  current  Jq,  and  the  other  Jj,  a  correction  to  Jq  due  to  all 
other  effects  not  included  in  PO. 

J  =  J0  *  Jj  (21) 

Ufimtsev  refers  to  these  two  currents  as  the  "uniform"  and 
"nonuniform"  components,  respectively.  Jq  is  the  GO  field  on  an 
infinite  flat  surface  tangent  to  the  body  at  the  point  in  question. 

The  Jj  is  due  to  perturbations  created  by  the  departure  of  the  surface 
from  an  infinite  plane.  In  his  monograph  [4]  and  other  publications, 
Ufimtsev  considers  only  the  nonuniform  component  due  to  edges  and  in 
practice,  therefore,  PTD  is  a  high  frequency  method  for  estimating 
scattering  by  edged  bodies. 

The  nonuniform  component  near  an  edge  is  approximated  by  the 
current  appropriate  to  an  infinite  wedge,  or  half-plane  tangent  to  the 
edge  and  the  fields  scattered  by  an  edged  structure  are  then  calculated 
by  integrating  both  types  of  current  over  the  surface. 

The  explicit  solution  for  the  nonuniform  component  is  extremely 
difficult  because  of  very  complicated  and  immense  equations  involved. 
The  contribution  of  the  nonuniform  component  is  determined  by  analogy 
with  the  corresponding  solution  for  a  strip.  This  leads  to  an 


expression  for  the  nonuniform  contribution  at  wide  angles  in  terras 
of  modified  diffraction  coefficients  representing  the  fields  scattered 
by  overlapping  half-planes  whose  orientation  is  specified  by  the  flash 
points  obtained  from  a  stationary  phase  analysis.  The  nonuniform 
portion  of  the  total  scattered  field  is  thereby  reduced  to  a  summation 
of  scattered  fields  rather  than  surface  currents. 


2.6  Numerical  Methods 

The  starting  point  for  the  numerical  method  is  the  vector 
integral  equation  for  the  surface  current.  It  is  the  most  logical 
starting  point  because  boundary  and  radiation  conditions  are  auto¬ 
matically  taken  into  account.  To  get  our  equation,  we  have  to 
rewrite  Eq.  (12a) 


y  J(r)  =  nx}?1  (r)  ♦  -j - 
r,  rl  e  S 


The  next  step  is  to  divide  the  scattering  surface  in 


portions  over  each  of  which  the  amplitude  and  the  phase  of  the  surface 
field  can  be  considered  as  approximately  constant  or  as  varying  in  an 


approximately  known  way.  Thus  the  integral  sign  in  Eq.  (22)  can  be 
replaced  by  a  summation  sign  and  we  get  a  set  of  linear  algebraic 
equations  to  be  solved  numerically  by  a  computer. 

The  set  of  algebraic  equations  can  be  reduced  to  a  matrix  form 


[V]  [J]  =  [H1] 


(23) 


where  [V]  is  the  known  coefficient  of  the  current  [J],  while 
[H1]  is  the  incident  magnetic  field.  So  [J]  can  be  found  out  by- 
inverting  matrix  [V]. 

Equation  (22)  can  be  solved  by  Moment  method.  For  details 
see  the  book  by  Harrington  [5],  A  short  review  will  be  given  here. 

The  idea  of  moment-method  solutions  is  that  the  linear  operator 
equation 

L(f)  =  g  (24) 

where  L  is  a  linear  operator,  g  is  some  known  driving  fn.  and  f  is 
to  be  determined,  may  be  solved  by  first  expanding  f  in  a  series  of 
basis  fns.  fR  with  unknown  coefficients  an>  viz, 

f  -  l  Vn  (“1 

n 

Next,  a  set  of  weighting  or  testing  fns.  U)m  is  chosen  and 
an  inner  product  taken  in  Eq.  (24).  After  substitution  of  Eq.  (25) 
into  Eq.  (24)  the  result  is 

Lfn>*<V«>  <26> 

which  may  be  written  as 

'W  t“„1  *  IsJ 

This  is  easier  to  invert,  and  matrix  inversion  yields  the 


solution  for  the  unknown  coefficients. 


After  finding  the  current  on  the  surface  we  can  use  Eq.  (16)  or 


(17)  to  find  the  scattered  field. 

The  numerical  method  is  very  effective  when  the  shape  of  the 
body  is  complex.  It  is  also  well  suited  to  problems  in  the  low 
frequency  and  resonance  regions . 

If  the  dimensions  of  the  scatterer  are  very  large  compared  to 
the  wavelength,  the  number  of  linear  algebraic  equations  become  very 
large  and  difficult  to  handle  by  presently  available  computers. 


in  which  ft  is  the  position  vector,  and  kg  is  the  phase  constant  of 
empty  space.  Substituting  (28)  into  (9)  and  (3)  and  equating  like 
powers  of  w,  one  obtains  the  eikonal  equation 


together  with  the  first-order  transport  and  conditional  equations 


3s  + 


1 

7 


&E 
n  o 


=  0 


a  •  t  =  0 

o 


(2 

(3 


plus  higher  order  transport  and  conditional  equations  which  do  not 
concern  us  here.  In  the  preceding  equations 

S  =  ^  ( 

n  v 

is  a  unit  vector  in  the  direction  of  the  ray  path,  and  s  is  the 
distance  along  the  ray  path. 

We  are  interested  here  in  the  solution  at  the  high  frequency 
limit,  so  the  asymptotic  approximation  of  E  reduces  to 

Its)  -  e-^^ijs)  (: 

Equation  (3TQ  is  readily  integrated  to  obtain 


£(S)  =  E  (0) 
0 


e-jk  *(0) 
e  o 


P  P 

_  1  2 

(p  +s) (p+s) 
]  2 


-jks 


in  which  s  =  0  is  taken  as  a  reference  point  on  the  ray  path,  and 
p  and  p  are  the  principal  radii  of  curvature  of  the  wavefront  at 


Employing  the  Maxwell  curl  equation 


V  x  £  =  -jkzft 


(7) 


we  can  obtain  H  . 

Equation  (34)  gives  the  propogation  of  £  in  the  space. 

Let  us  find  the  GO  field  reflected  from  the  point  QR  on  a 
perfectly-conducting  smooth  curved  surface  S;  the  distance  between 
Qr  and  the  field  point  on  the  reflected  ray  is  denoted  by  s.  The 
outward  directed  unit  normal  vector  at  QR  is  n,  and  s  and  s  are 
unit  vectors  in  the  directions  of  incidence  and  reflection,  respec¬ 
tively,  as  shown  in  Fig.  2. 

From  the  boundary  condition  on  the  total  electric  field  at 

Qr  on  S 


ft  x  (F+F)  =  0 


(35) 


it  can  be  shown  that 


^(Qo) 


;  ■  ?>r> 


[8 


i  g  r 
II  II 


-  e. 


i' 


(36) 


where  e^ 
and  , 
so  that 


is  a  unit  vector  perpendcular  to  the  plane  of  incidence, 
e^1"  are  unit  vectors  parallel  to  the  plane  of  incidence 


x  s 


(37) 
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In  matrix  notation  the  reflection  coefficient  has  a  form 


familiar  for  the  reflection  of  a  plane  electromagnetic  wave  from  a 
plane,  perfectly  conducting  surface,  namely 


R  = 


(3 


This  is  not  surprising  if  one  considers  the  local  nature  of  high- 
frequency  reflection,  i.e.,  the  phenomenon  for  the  most  part  depends 
on  the  geometry  of  the  problem  in  the  immediate  neighborhood  of  QR. 
Note  that  the  incident  and  reflected  fields  must  be  phase  matched  on 
S  to  satisfy  Eq.  (35),  i.e., 

/(Qr)  »  *r(QR)  (3! 


The  above  equality  leads  to  the  law  of  reflection,  and  it  is  also 
employed  to  obtain  (36) . 

The  geometrical  optics  reflected  ray  is 

=  £lg(QR)’P- 


4 


-jks 


(pr+s)  (Pr+s) 
1  2 


(41 


p  r  and  pr  are  the  principal  radii  of  curvature  of  the  reflected 
1  2 

wavefront  at  QR.  The  general  expressions  have  been  given  by 
Kouyoumj ian  [6]. 

For  our  study  we  are  going  to  assume  that  the  incident  wave 

T  T 

is  a  plane  wave,  and  s  >>  p  ,  p  so  that 

1  2 


pT  pr 


a  a  /4 


(4 


and  Eq.  (40)  becomes 


F«cs)  =  f1* (qr)-r  /rr  e-Jks 
R  -  *  2  ~rr 


a  and  a  are  principal  radii  of  curvature  of  the  surface  S  at  Qn. 
1  2  R 


It  should  be  noted  that  for  u>  sufficiently  large  the  GO  field 


requires  no  correction,  i.e.,  the  scattering  phenomenon  is  entirely 


dominated  by  GO,  This  is  the  case  for  backscatter  from  smooth  curved 


surfaces  with  radii  of  curvature  very  large  in  terms  of  a  wave  length 


2.8  Geometrical  Theory  of  Diffraction  (GTD) 


Geometrical  Optics  fails  to  predict  the  scattering  from  edges 


and  also  in  the  shadow  region.  To  overcome  these  limitations.  Dr. 


Keller  introduced  a  theory  which  is  now  known  as  Geometrical  Theory 


of  Diffraction.  Keller  [7-9]  has  shown  how  the  diffracted  field  may 


be  included  in  the  high  frequency  solution  as  an  extension  of  GO.  The 


postulates  of  GTD  are  summarized  as  follows. 


(1)  The  diffracted  field  propagates  along  rays  which  are 


determined  by  a  generalization  of  Fermat's  principle  to  include  points 


on  the  boundary  surface  in  the  ray  trajectory. 


(2)  Diffraction  like  reflection  and  transmission  is  a  local 


phenomenon  at  high  frequencies,  i.e.,  it  depends  only  on  the  nature 


of  the  boundary  surface  and  the  incident  field  in  the  immediate  neigh¬ 


borhood  of  the  point  of  diffraction. 


(3)  The  diffracted  ray  propagates  along  its  ray  such  that 
(a)  power  is  conserved  in  a  tube  (or  strip  of  rays) 


(b)  The  phase  delay  along  the  ray  path  equals  the  product 
of  wave  number  of  the  medium  and  the  distance. 

Diffracted  rays  are  initiated  at  points  on  the  boundary  surface 
where  the  incident  GO  field  is  discontinuous.  Examples  of  such  points 
are  edges,  vertices  and  points  at  which  the  incident  ray  is  tangent  to 
a  smooth  curved  surface.  The  diffracted  rays  like  the  GO  rays  follow 
paths  which  make  the  optical  distance  between  the  sourcepoint  and 
the  field  point  an  extremum,  usually  a  minimum. 

Let  us  consider  the  field  radiated  from  a  point  source  at  0, 
and  observed  at  P  in  the  presence  of  a  perfectly  conducting  wedge, 
as  shown  in  Fig.  3,  where  the  rays  are  projected  on  a  plane  perpendic¬ 
ular  to  the  edge  at  the  point  of  diffraction  Q£. 

Applying  the  generalized  Fermats  principle,  the  distance  along 
the  ray  path  OQgP  is  a  minimum  and  the  law  of  edge  diffraction 

s  •  e  =  s  •  e  (43) 

results.  Here  §  is  a  unit  vector  directed  along  the  edge,  and  §'  and 
s  are  unit  vectors  in  the  directions  of  incidence  and  diffraction, 
respectively.  Eq.  (43)  also  follows  from  the  requirement  that  the 
incident  and  diffracted  fields  be  phase  matched  along  the  edge.  If 

the  incident  ray  strikes  the  edge  obliquely,  making  an  angle  &' 

/  ° 
with  the  edge,  as  shown  in  Fig.  4,  diffracted  rays  lie  on  the  surface 

of  a  cone  whose  half  angle  is  equal  to  The  position  of  the 

diffracted  ray  on  this  conical  surface  is  given  by  angle  <f>,  and  the 

direction  of  the  ray  incident  on  the  edge,  by  the  angles  <p'  and 


It  has  been  shown  by  Kouyoumjian  [6]  that  the  expression  for 


I 

} 

*  the  electric  field  of  the  edge  diffracted  ray  is 

J'w  ■  e;  )  J  e-jfcs  (44) 

in  which  D (4> , 4> ^ , 6 " D  is  the  dyadic  edge  diffraction  coefficient,  p  is 
-  o 

f  the  distance  between  the  caustics,  one  is  the  point  of  diffraction 

itself  and  the  second  caustic  of  the  diffracted  ray  is  away  from  the 
surface.  The  dyadic  diffraction  coefficient  for  a  perfectly  conduct- 
|  ing  wedge  has  been  obtained  by  Kouyoumjian  and  Pathak  [10].  They 

found  out  D  from  the  asymptotic  solution  of  canonical  problems,  which 

♦ 

in  this  case  involves  the  illumination  of  the  wedge  by  plane,  cylind- 
,  rical,  conical,  and  spherical  waves.  The  solution  of  these  canonical 

r 

|  problems  serves  as  a  basis  for  deducing  the  dyadic  diffraction  coeffi- 

|  cient  for  arbitrary  wavefront  illumination  and  for  the  more  general 

case  where  there  are  curved  edges  and  curved  surfaces. 

* 

i 

Let  us  introduce  an  edge-fixed  plane  of  incidence  containing 
the  incident  ray  and  the  edge  and  a  plane  of  diffraction  containing 
the  diffracted  ray  and  the  edge.  The  unit  vectors  and  $  are 
perpendicular  to  the  edge  fixed  plane  of  incidence  and  the  plane  of 

A  A 

diffraction,  respectively.  The  unit  vectors  6  '  are  8  parallel  to 

o  o 

the  edge  fixed  plane  of  incidence  and  the  plane  of  diffraction 
respectively,  and 

8  3'*  ,  8  =  s  x  (45) 

0  ’9 
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Thus  the  coordinates  of  the  diffracted  ray  (s,B  ,<p )  are  spherical 

0 

coordinates  and  so  are  the  coordinates  of  the  incident  ray 

(s',  8  ' ,  4>')  except  that  the  incident  (radial)  unit  vector  points 
o 

toward  the  origin  Q^.  These  ray  fixed-fixed  coordinates  and  their 
unit  vectors  are  shown  on  Fig.  5  and  Fig.  6. 

If  the  dyadic  diffraction  coefficient  is  expressed  in  an  edge 
fixed  coordinate  system,  it  is  found  to  be  the  sum  of  seven  dyads.  In 
matrix  notation  this  means  that  the  diffraction  coefficient  is  a 
3x3  matrix  with  seven  non-vanishing  elements,  instead  of  the  2x2 
diagonal  matrix  which  may  be  used  to  represent  the  diffraction  coeffi¬ 
cient  in  the  ray-fixed  coordinate  system.  In  this  sense  the  ray- 
fixed  coordinate  system  is  the  natural  coordinate  system  of  the  problem. 

For  each  type  of  illumination  mentioned  previously,  it  is  shown 
by  Pathak  and  Kouyoumjian  [11]  that  the  dyadic  diffraction  coefficient 
can  be  represented  simply  as  the  sum  of  two  dyads,  if  the  ray-fixed 
coordinates  mentioned  in  the  preceding  paragraph  are  used. 

D(<J>,4>',8  ')  =-8'  8  D  f4>,4>',6')  -  OD h(<M', S')  (46) 

where  Ds  is  the  scalar  diffraction  coefficient  for  the  acoustically 

soft  (Dirichlet)  boundary  condition  at  the  surface  of  the  wedge,  and 

D,  is  the  scalar  diffraction  coefficient  for  the  acoustically  hard 
h 

(Neumann)  boundary  condition. 

Expressions  for  the  scalar  diffraction  coefficients  which  are 
valid  at  all  points  away  from  the  edge  are 


2n ./  2-rrk 
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\  cot  (— 2^j— )  F (kL  a  (a  )) 


Ds  " 
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sin  6 


♦  r 

♦  cot  (!£-)  F(kLa"(a"))  +  [cot  (—-)  F(U  V(a+)) 


♦  cot  (2^2_)  FfU^Vfo*))]} 


where  a  =  <p  + 


a  -  <p  -  (p' 


F (X)  =2j|/Tl  e"^X 


00 

'  2 


e"jT  dx 


|/Tl 


(47) 

(48) 

(49) 

(50) 


If  X  is  greater  than  10  then  F(X)  -  1.0. 

.  +  + 

The  large  argument  of  the  transition  fq.X  *  kL  a~(a  ) 
is  KL^T°  ,  where  L^10  is  a  distance  parameter  which  will  be  given 


later. 


Let  <p 


=  a 


then 


— ,  — s  -  (.2n'nN  -  ct  , 

a  (a  )=  2  cos  ( - = - ) 


(51) 


* 

in  which  N  are  the  integers  which  most  nearly  satisfy  the  equations 


+ 


2imN*  -  a  =  n  (52) 

27mN"  -  a  *  -it  (53) 


a  (a  )  is  a  measure  of  the  angular  separation  between  the 
field  point  and  a  shadow  or  reflection  boundary.  The  +  or  -  super¬ 
scripts  are  ssociated  with  the  integers  N+  and  N  ,  respectively, 
which  are  defined  by  Eq.  (52)  and  (53).  For  exterior  edge  diffraction 
(l<n<2)  N*  =  0,  or  1,  and  N~  =  -1,  0,  or  1. 

The  distance  parameter  L  can  be  found  by  imposing  the  condition 
that  the  total  field  (the  sum  of  GO  field  and  the  diffracted  field) 

be  continuous  at  shadow  or  reflection  boundaries. 

,  •  .  * 
p1+Sjp1  p1  sin28 


P.(P  +s)(p  +s) 
e  i  2 


,pr+s>  pr  pr  sin2  6 


P  (P  +s)  (P  +s) 

1  2 

wher  pr  and  pr  are  the  principal  radii  of  curvature  of  the  reflected 
1  2 

■»» 

wavefront  at  Q^,  and  p  is  the  distance  between  the  caustics  of  the 

diffracted  ray  in  the  direction  of  the  reflection,  p1  and  p1  are  the 

12 

principal  radii  of  curvature  of  the  incident  wavefront  at  Qg,  and 

p  1  is  the  radius  of  curvature  of  incident  wavefront  in  the  edge  fixed 
e 

plane  of  incidence.  The  additional  superscripts  0  and  n  on  L  in  Eq. 
(47)  denote  that  the  radii  of  curvature  (and  caustic  distance  p)  are 
calculated  at  the  refelection  boundaries  and  (2n- 1 ) tt- "  ,  respec¬ 

tively.  In  the  far  zone  where  s  »  p  and  the  principal  radii  of  curva¬ 
ture  p  and  p  of  the  incident  and  reflected  wavefronts  at  QP  Eq.  (54) 
2  2  c» 


and  (55)  simplify  to 


p  p  sin2  Bn” 


2 


The  appropriate  superscripts  are  omitted  here  for  the  sake 
of  simplicity. 

The  cuastic  distance  p  can  be  calculated  employing  differential 
geometry.  It  is 


1  1 
p  =  ^r 


fie*(s-S) 
a  sin2S 


(57) 


ft  is  the  unit  vector  normal  to  the  edge  at  QP  and  directed  away  from 
the  center  of  curvature,  a  is  the  radius  of  curvature  of  the  edge  at 
<JE,  a  >  0. 

In  summary,  to  calculate  the  diffraction  from  a  curved  edge, 
the  scalar  diffraction  coefficients  from  Eq.  (47)  are  substituted 
into  Eq.  (46)  and  the  resulting  dyadic  diffraction  coefficient  is 
substituted  into  Eq.  (44).  The  cuastic  distance  p  is  calculated  from 
Eq.  (57).  In  matrix  notation  the  diffracted  electric  field  takes  the 
form 


(58) 


where  D(0,  <$>,  O',  <p')  is  the  vertex  diffraction  coefficient.  Very  little 
work  has  been  done  on  the  high-frequency  diffraction  by  vertices. 


When  an  incident  ray  strikes  a  smooth,  curved  perfectly  conduct¬ 
ing  surface  at  grazing  incidence,  i.e.,  at  the  shadow  boundary,  a 
part  of  the  energy  is  diffracted  into  the  shadow  region.  Let  us 
consider  the  field  radiated  by  the  source  0  and  observed  at  p  in  the 
shadow  region  as  shown  in  Fig.  7  and  Fig.  8. 

To  satisfy  the  generalized  Fermats  principle,  the  surface 
ray  must  be  a  geodesic  curve.  Also  it  is  required  that  the  incident 
and  diffracted  fields  are  phase  matched  to  the  surface  ray  field  at 
Qj  and  C^.  It  is  assumed  that  the  surface  rays  are  torsionless,  i.e., 
£  does  not  change  direction  along  the  surface  ray.  Kouyoumjian  [6] 
has  given  the  following  expression  for  surface  diffraction. 

-  ^CQ,)  -  Ifi.V  *  bU]  yj-^  e'jl<s  (60) 

In  the  same  reference  Kouyoumjian  [6]  has  shown  how  to  find  F 
and  G  and  p  is  found  from  Eq.  (57).  Calculating  F  and  G  is  quite 
involved,  and  interested  readers  are  requested  to  see  reference  [6]. 
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SIDE  VIEW 


CHAPTER  THREE 


THE  CHOICE  OF  METHOD 

Recently  there  has  been  considerable  interest  in  the  integral 
equation  formulation  of  the  radiation  problem,  and  its  solution  by 
the  moment  method.  Arbitrary  shapes  can  be  handled  by  this  method, 
but  in  general,  numerical  results  are  restricted  to  objects  not  large 
in  terms  of  a  wavelength,  because  of  the  limitations  of  present  day 
computers. 

Tsai,  Dudley  and  Wilton  [12]  have  used  the  moment  method  to 
solve  the  vector  magnetic  field  integral  equation,  and  thus  calculated 
electromagnetic  scattering  by  a  three  dimensional  rectangular  box. 

They  have  shown  that  using  an  8x8  sampling  on  each  face  of  a  cube 
would  require  768x768  matrix  to  be  solved.  Looking  at  this  enormous 
problem  for  such  a  simple  object  as  a  cube,  we  shelve  this  method  for 
the  time  being  until  some  mathematical  breakthrough  or  the  computer 
improvement  is  achieved  in  solving  such  problems. 

Now  we  are  left  with  two  methods  for  calculating  the  scattered 
field.  One  is  Physical  Optics  combined  with  Physical  Theory  of 
Diffraction  and  the  other  is  Geometrical  Optics  combined  with  Geo- 
meticral  Theory  of  Diffraction.  The  results  obtained  by  PO  are 
frequency  dependent,  and  it  is  therefore  possible  that  PO  provides  a 
more  accurate  estimate  of  the  scattering.  But  it  has  been  noted  that 
PO  in  general  fails  to  satisfy  the  reciprocity  theorom  everywhere 
except  in  the  direction  of  a  specular  return.  Also  no  rigorous  justi¬ 
fication  for  PO  approximations  have  been  established.  PTD  has  a  lot 


of  potential  to  be  a  leading  method  in  predicting  diffraction,  but 
a  lot  of  work  will  have  to  be  done.  Thus  if  PTD  is  to  be  judged  by 
the  description  given  by  Ufimtsev  [4],  it  is  less  complete  than  GTD 
and,  moreover,  does  not  give  the  correct  second  order  diffraction  for 
strips  and  disks. 

According  to  Knott,  Thomas  and  Senior  [13]  GTD  is  as  good  as 
any  other  method.  It  has  some  advantages  and  some  disadvantages. 

The  serious  disadvantage  being  that  it  predicts  infinite  field  at  the 
ray  caustic.  So  as  long  as  we  keep  away  from  the  ray  caustic,  it  is 
an  elegant  method  with  the  conceptual  simplicity.  Also  R.  Mittra 
and  W.  L.  Ko  [14]  have  concluded  in  this  reference  that  GO  solution 
is  a  very  good  approximation  to  the  exact  solution  except  in  the 
neighborhood  of  the  forward  scattering  direction.  Also,  it  has  been 
found  in  many  cases  GO  and  GTD  has  worked  surprisingly  well  on 
radiating  objects  as  small  as  a  wavelength  or  so  in  extent. 

Since  we  are  interested  only  in  the  backscattering,  we  conclude 
that  at  the  present  time  the  best  practical  method  to  model  the  scatt¬ 
ering  field  from  a  3-D  object  is  the  Geometrical  Optics  in  combination 
with  Geometrical  Theory  of  Diffraction. 


CHAPTER  FOUR 
APPLICATIONS 

Application  of  Geometrical  Optics  and  the  Geometrical  Theory 
of  Diffraction  to  calculate  the  backscattered  field  from  3-D  objects 
is  a  straightforward  matter.  The  first  step  is  to  find  out  all 
scattering  points  on  the  object.  This  is  purely  an  exercise  in 
geometry. 

Since  we  are  interested  only  in  the  backscattered  field,  it  is 
easier  to  find  the  scattering  points.  The  direction  of  the  scattered 
ray  must  be  exactly  opposite  to  the  incident  ray.  The  points  which 
satisfy  this  constraint  are  the  specular  points. 

The  second  step  is  to  apply  the  necessary  equations  that  are 
developed  before,  at  these  scattering  points.  The  final  step  is  to 
sum  up  all  the  scattered  field  at  the  given  position  to  get  the 
total  field. 

We  are  going  to  apply  this  methodology  to  calculate  the  back- 
scattered  field  from  a  sphere  which  has  been  studied  extensively  and 
from  a  spherical  segment. 

4.1  Back  Scattering  from  a  Sphere 
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A  great  deal  of  research  has  been  done  on  the  EM  backscatter 
from  the  sphere  since  it  has  the  simplest  geometry.  The  exact  closed 
form  mathematical  solution  for  the  EM  backscatter  is  known.  A  theor¬ 
etical  derivation  of  the  field  scattered  by  perfectly  conducting 
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spheres  is  given  by  Harrington  [15].  Chapter  10  of  Electromagnetic 
and  Acoustic  Scattering  by  Simple  Objects  [16]  gives  the  pertinent 
equations  to  calculate  the  exact  EM  backscatter  from  spheres. 
Rheinstein  [17]  and  Bechtel  [18]  have  tabulated  the  numerical  results. 

The  purpose  of  picking  up  the  sphere  to  calculate  the  back- 
scattering  is  to  test  the  practicality  of  our  equations  with  the  known 
results.  The  situation  is  shown  in  Fig.  9. 

In  this  case  we  will  have  backscattering  due  to  reflection  and 
due  to  surface  diffraction.  But  we  will  neglect  the  contribution  of 
the  surface  diffraction  since  it  is  quite  small  compared  to  that  of 
reflection.  Dimensional  considerations  show  that  the  coefficient  of 
surface  diffraction  decreases  exponentially  with  decreasing  X,  and 
our  mathematical  model  hinges  on  the  assumption  that  the  frequency  of 
EM  radiation  is  very  high.  Therefore,  at  high  frequencies  the  return 
from  the  surface  wave  will  be  very  weak. 

The  specular  point  will  lie  at  the  intersection  of 
the  surface  of  the  sphere  and  the  line  joining  the  center  of  the 
sphere  and  the  aircraft  radar  for  the  obvious  reason  that  the  incident 
ray  will  be  orthogonal  to  the  surface,  so  that  the  reflected  ray  will 
be  exactly  in  the  opposite  direction. 

The  equation  to  be  applied  is  the  GO  formulation  and  we  get 

=  4 

0 

F«(s)  is  the  scattered  field  received  at  the  antenna,  i=  is  the  EM 
field  af*the  origin  of  antenna. 


p  and  p  are  the  principal  radii  of  curvature  of  the  wavefront  at 
12 

the  origin  of  the  antenna  and  it  is  assumed  to  be  known,  a^  and  a^ 
are  the  principal  radii  of  curvature  of  the  surface  at  the  point  of 
reflection,  and  in  our  case  it  is  equal  to  the  radius  of  the  sphere, 
s  is  the  distance  between  the  antenna  dnd  the  specular  point,  k  is 
the  wave  number. 

4.2  Back  Scattering  from  a  Spherical  Segment 

Now  we  are  going  to  calculate  the  backscattered  field  using 
GO  and  GTD , 

)  The  situation  that  we  will  be  considering  is  as  follows.  An 

airplane  in  the  sky  is  trying  to  detect  objects  like  cars,  tanks, 
etc.  which  are  on  the  ground  with  its  radar.  The  radar  emits  the  EM 
signals  and  then  gathers  the  scattered  signals.  Since  the  velocity 
of  the  plane  is  much  lower  than  that  of  EM  signals,  we  can  assume 
that  the  whole  situation  is  static  and  the  signals  that  the  radar 
receives  are  the  backscattered  signals.  We  will  also  assume  the 
radar  to  be  a  point  source  and  a  point  receiver  because  the  distance 
between  the  object  and  the  radar  is  much  larger  than  the  dimensions 
of  the  radar.  We  will  also  assume  that  the  object  is  perfectly 
conducting.  Let  us  take  the  object  to  be  a  car.  To  make  the  problem 
tractable,  we  will  say  the  car  resembles  a  spherical  segment  sitting 
on  the  ground.  The  situation  has  been  depicted  in  Fig.  10. 

To  summarize,  the  simplifications  that  we  have  made 
(i)  The  object  is  perfectly  conducting. 
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(ii)  The  object  is  a  spherical  segment. 

(iii)  The  object  is  very  large  compared  to  the  wavelength  of 
the  radiation. 

(iv)  The  radiation  is  mono chromatic  and  the  surrounding 
of  the  object  is  almost  a  vacuum. 

(v)  The  frequency  of  the  signal  is  quite  high. 

(vi)  The  situation  is  static. 

We  will  not  consider  backscattering  due  to  surface  diffraction 
for  the  following  reason.  Dimensional  considereations  show  that  the 

j_ 

2 

edge  diffraction  coefficients  are  proportional  to  X  and  tip  or  vertex 
diffraction  coefficients  to  X.  The  field  diffracted  around  a  curved 
surface  decreases  exponentially  with  decreasingX  and  is  consequently 
weaker  than  that  diffracted  by  an  edge  for  low  values  of  X  [9] .  This 
is  even  more  true  in  the  case  of  backscattering.  To  substantiate  our 
claim  further,  we  refer  to  the  research  done  by  Weiner  [19].  He  con¬ 
fined  his  studies  only  to  ahemi sphere  with  a  nose  on  incidence.  He 
also  conducted  some  experiments.  He  concluded  that  the  GTD  appeared 
to  give  the  correct  magnitude  and  frequency  dependence  for  the  singly 
diffracted  return.  However,  the  doubly  diffracted  return  is  consider¬ 
ably  weaker  than  diffraction  theory  would  predict  and  decreases 
rapidly  with  frequency.  Surface  wave  return  is  the  weakest  of  all  and 
it  also  decreases  rapidly  with  the  increasing  frequency. 

The  origin  of  the  coordinate  axis  coincides  with  the  center  of 
an  imaginary  sphere  of  which  the  spherical  segment  is  a  part.  Let  D 
be  the  diameter  of  the  imaginary  sphere  and  let  the  height  of  the 


spherical  segment  be  h.  The  flat  surface  of  the  spherical  segment 
lies  parallel  to  the  XY  plane. 


When  the  radiated  beam  covers  the  object,  we  will  have  scatter¬ 
ing  from  the  object.  At  all  times  there  will  be  at  least  one  scatter¬ 
ing  point  on  the  edge.  Additional  scattering  points  will  come  into 
the  picture  depending  upon  the  relative  position  of  the  aircraft  and 
the  spherical  segment.  In  one  situation  (see  Fig.  12)  the  scattering 
will  be  from  the  edge  only. 

In  another  situation  (see  Fig.  12)  there  will  be  two  scattering 
points — one  point  will  be  on  the  edge,  while  the  other  is  on  the 
smooth  curved  surface. 

In  the  third  situation  (see  Fig.  13)  there  will  be  three 
scattering  points,  two  points  are  located  on  the  edge  and  the  third 
point  is  on  the  smooth  curved  surface.  We  will  not  consider  the  nose 
on  incidence  because  it  is  a  remote  case. 

First  let  us  give  the  conditions  under  which  these  additional 
scattering  points  occur. 

Let  the  positional  vector  of  the  aircraft  be  it 

$  =  Xi  ♦  Wj  +  Hk  (62) 

Therefore,  the  angle  a  between  the  line  joining  the  aircraft 
with  the  origin  and  its  projection  on  the  gound  (XY  plane)  is  given 
by 
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The  angle  8  between  the  line  joining  any  point  on  the  curved 


edge  with  the  origin  and  its  projection  on  the  ground  is  given  by 

8  =  Arc  sin  0~-  j  (64) 

Whenever  the  radiation  covers  the  object  we  will  have  back- 
scatter  at  least  from  a  point  on  the  curved  edge  which  is  closest  to 
the  aircraft.  This  is  so  because  the  incident  radiation  must  be 
perpendicular  to  the  edge  in  order  to  get  the  backscatter. 

We  will  have  reflection  from  the  curved  surface  of  the  spheri¬ 
cal  segment  in  addition  to  the  diffraction  from  the  edge  when 

8  <  a  <  tt/2  (65) 

The  point  of  reflection  will  lie  at  the  intersection  of  the 
curved  surface  of  the  spherical  segment  and  the  line  joining  the 
radar  with  the  origin;  i.e.,  the  incident  and  the  reflected  ray  are 
perpendicular  to  the  surface. 

We  can  have  one  more  scattering  point  on  the  curved  edge.  This 
scattering  point  must  lie  diametrically  opposite  to  the  first  scatter¬ 
ing  point  because  the  incident  ray  must  be  orthogonal  to  the  tangent 
of  the  ecge  at  that  point  to  get  the  backscatter. 

Let  X  be  the  angle  between  the  line  joining  the  aircraft  with 
point  A  and  its  projection  on  the  ground. 

/  H-§*H  \ 

Y  *  arc  tanl  — ^ -  I  (66) 


where  Vi'  = 


(67) 


sj  X2+W2  ’♦  a' 

a"  is  the  radius  of  curvature  of  the  edge  and  it  is  given  by 

a'  =  v/ (D/2) 2  - (D/2-h) 2  (68) 

We  will  have  this  additional  scattering  point  when 


(J  -8)  <  Y 


(69) 


This  can  clearly  be  seen  from  Fig.  13. 

We  conclude  that  if  the  condition  given  by  Eq.  (69)  is  satis¬ 
fied,  we  will  have  diffraction  from  two  points  on  the  curved  edge 
and  reflection  from  one  point  on  the  surface.  If  only  Eq.  (65)  is 
satisfied,  then  we  will  have  scattering  from  one  point  on  the  edge 
and  the  reflection  from  the  surface.  Note  that  if  Eq.  (69)  is 
satisfied,  Eq.  (65)  is  satisfied  automatically.  If  the  condition  in 
Eq.  (65)  is  not  satisfied,  then  the  scattering  will  be  from  one  point 
on  the  edge  only. 

Let  us  consider  the  situation  where  there  is  diffraction  from 
two  points  on  the  edge  and  the  reflection  from  one  point.  Refer  to 
Fig.  14. 

We  have  the  solution  for  the  reflection  from  point  B 
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(70) 
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t®(0)  is  the  electric  field  returned  to  the  aircraft  antenna  after 
B 

the  reflection  at  point  B. 
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E  is  the  originating  field  at  the  radar. 

p  and  p  are  the  principal  radii  of  curvature  at  the  point  of  origin 
1  2 

of  radiation  and  they  are  assumed  to  be  known. 

s_  is  the  distance  between  the  antenna  and  the  point  B  . 

D 


a  and  a  are  the  radii  of  curvature  of  the  reflected  wavefront  at  B. 
1  2 


a  =  a  =  D/2 
1  2 


(72) 


k  =  2u/X 


(73) 


where  A  is  the  wavelength. 

Now  let  us  consider  the  diffraction  from  point  C  which  is  the 
closest  point  on  the  edge  from  the  antenna.  Let  ^  be  the  positional 
vector  of  point  C  with  respect  to  the  origin. 

Rc  =  j  (cosB  (cos0  i  +  sin0  j)  +  sing  k]  (74) 

0  =  arc  tan  (^)  (75) 

Consider  Fig.  15  and  Fig.  16. 

Let  us  introduce  an  edge  fixed  plane  of  incidence  containing 
the  incident  ray  and  the  tangent  of  the  edge  at  point  C.  Also  con¬ 
sider  the  edge  fixed  plane  of  diffraction  containing  the  diffracted 
ray  and  the  tangent  of  the  edge.  Since  the  incident  ray  and  the 
diffracted  ray  lie  on  the  same  line,  both  the  planes  are  identical. 
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in  > 


The  unit  vectors  <j>  and  <$>  are  perpendicular  to  the  edge  fixed  plane. 

A 

is  associated  with  the  incident  ray  while  4>  is  associated  with  the 
diffracted  ray. 


A  A  A  A  A 

4>  =  <t>'  =  sin(4>-ir)  [-cos0  i  -  sinQ  j]  +  cos(<j)-n)k  (76) 


+  7t  (77) 


A  A 

The  unit  vectors  8  '  and  8  are  parallel  to  the  edge  fixed  plane  of 
o  o 

incidence  and  diffraction.  Also  they  are  orthogonal  to  the  incident 

A 

and  diffracted  ray.  e  is  the  unit  vector  tangential  to  the  curved 
edge  at  the  point  C. 
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A  A 

-sin0  i  ♦  cos0  j 


(78) 

(79) 

(80) 


is  the  unit  directional  vector  of  the  diffracted  ray. 

The  angle  between  the  two  surfaces  which  form  the  edge  is 


given  by 

nTr  =  2tt-  (  j  -  Arc  Sin  (~07.h)  j 

n  _  3  1  Arc  Sin  (^p-) 

n  -  7  4  ?  0/2 


(81) 

(82) 


(C)  =  The  component  of  the  electric  field  vector  at  point  C 

A 

in  the  direction  of  B' 
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s^,  is  the  distance  between  the  point  C  and  the  antenna. 
sc  =  yj [X- j  cosB  cos0]2  +  [W-y  cosB  sin©]2  +  [H-R  sinB]2 

E*^(C)  =  The  component  of  the  incident  electric  field  vector  at 

<P 

point  C  in  the  direction  of  <J> 

(OC)  =  The  component  of  the  electric  field,  diffracted  from 

g  A 

point  C,  returned  to  the  antenna  in  Bfl  direction. 
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where  Dg  is  called  the  diffraction  coefficient  for  the  soft  surface 
and  p  is  caustic  distance. 


p  is  given  by 
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p  1  is  the  radius  of  curvature  of  the  incident  wavefront  in 
e 

the  edge  fixed  plane  of  incidence.  Since  the  wavefronts  are  almost 
planar  at  point  C 


=  0.0 


n  is  the  unit  vector  normal  to  the  edge  at  C  and  directed 
e 

away  from  the  center  of  curvature. 


Therefore 


(s'-s)  =  -2  cos 


The  component  of  the  electric  field,  diffracted 

a 

from  point  C,  returned  to  the  antenna  in  4>  direction. 
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(90) 


is  called  the  diffraction  coefficient  for  the  hard 


surface. 

The  general  expression  for  Ds  and  is  given  by  Eq.  (47). 
In  this  particular  case 
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(91) 


The  upper  sign  goes  with  Ds,  while  the  lower  sign  goes  with  D^. 
F(*)  is  called  a  transition  fn.  and  for  our  purpose 

F (• )  =  1.0  (92) 


The  transition  fn.  is  important  at  only  very  few  points, 
especially  when 


tt+2$ 

2n 


IT 


(93) 


At  this  particular  value  F(*)  has  different  value  than  1.0,  but 
this  is  a  solitary  point.  Therefore,  for  most  practical  case 
F(-)  =  1.0  . 
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Then  Eq.  (91)  simplifies  to 


-  <  -tt+24)..  .  ,tt-2<J>,  \ 

*  i cot  (-sr>  * cot  '-ir5 } 

(94) 

For  the  point  A  exactly  the  same  procedure  is  to  be  used, 
except  for  the  minor  changes  in  the  angles  and  distances. 

Let  be  the  positional  vector  of  point  A  which  is  the 
farthest  point  on  the  edge  from  the  antenna. 

-v  D  A  A  ^ 

Ra  =  7  [cos8(-cos9  i  -sin0  j)  +  sinB  k]  (95) 

A  A 

Let  <j>'  and  4>A  be  the  unit  vectors  perpendicular  to  the  edge 

A 

fixed  plane  at  point  A.  <J>^  is  associated  with  the  incident  ray  at  A, 
while  0A  is  associated  with  the  diffracted  ray. 

S'  =  *  (96) 

This  is  so  because  the  radius  of  curvature  of  the  edge  is 

much  smaller  than  the  distance  between  the  object  and  the  antenna. 

Also  since  the  incident  ray  and  the  diffracted  ray  have  exactly 

opposite  direction,  the  edge  fixed  planes  are  the  same. 

The  unit  vectors  §'  and  §  are  parallel  to  the  edge  fixed 

o  o 

plane  of  incidence  and  diffraction  and  they  are  the  same  as  before. 

£g'(A)  =  The  component  of  incident  electric  field  vector  at 
o  ~ 

point  A  in  the  direction  of  0' 
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-H  +  -  /  P  P  -  jks.  A 

Ej.(A)  =  (E  .0')  \/  7- - ^ - *  e  A  6' 
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is  the  distance  between  the  point  A  and  the  aircraft  antenna. 
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^-(A)  =  The  component  of  incideht  electric  field  vector  at 
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point  A  in  the  direction  of 
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z&  (°A^  =  Component  the  diffracted  electric  field  vector 


at  the  antenna  in  0  direction. 
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PA  is  the  distance  between  the  caustic 


PA  =  2  cos  (4>  >_tt) 


(101) 


<p'=  2 it  -  Arc  tan 
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=  The  component  of  the  electric  field  at  the  antenna, 
after  the  diffraction  at  point  A,  in  the  direction  of  4>  . 
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The  edge  diffraction  cbefficient  D'  and  D'  is  given  by 


-e 


-jkn/4 


^  2n  y/  2irk 


2  sinn/n  + 

cosir/n-l 


.tr+24)'.  .  ,71-24)  \ 

cot^-2ni  +  cot(~l*-> 


(104) 


The  total  scattered  field  returned  to  the  antenna  is  given  by 
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CHAPTER  FIVE 


HIERARCHY  PLUS  INPUT- PROCESS-OUTPUT  (HIPO)  CHARTS 

We  need  to  document  the  method  that  has  been  developed  in 
order  to  facilitate  an  overview  of  the  system  and  to  determine  program 
functions  for  coding  purposes.  One  such  technique  used  is  called 
Hierarchy  plus  Input- Process -Output  (HIPO)  which  is  developed  by  IBM. 

A  typical  HIPO  package  contains  three  kinds  of  diagrams. 

They  are 

(i)  Visual  Table  of  Contents:  This  diagram  contains  the 
names  and  identification  numbers  of  all  the  overview  and  detail  HIPO 
diagrams  in  the  package  and  relationship  of  the  functions  in  a 
hierarchical  fashion. 

(ii)  Overview  diagrams:  High  level  HIPO  diagrams,  called 
overview  diagrams,  describe  the  major  functions  and  reference  the 
detail  diagrams  needed  to  expand  the  functions  to  sufficient  detail. 

The  overview  diagrams  provide,  in  general  terms,  the  inputs,  processes, 
and  outputs. 

(iii)  Detail  diagrams:  Lower  level  HIPO  diagrams  contain  the 
fundamental  elements  of  the  package.  They  describe  the  specific 
functions,  show  specific  input  and  output  items,  and  refer  to  other 
detail  diagrams. 

For  more  detail  information  about  HIPO  charts,  please  refer  to 
IBM  manual  [20] . 

See  HIPO  package  in  Appendix  I. 
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CHAPTER  SIX 


DISCUSSION  AND  SUGGESTIONS 

We  discussed  the  merits  and  shortcomings  of  the  three  more 
popular  methods  of  calculating  scattering,  namely  Physical  Optics  (PO) 
together  with  Physical  Theory  of  Diffraction  (PTD) ,  Numerical  Methods, 
and  Geometrical  Optics  (GO)  in  conjunction  with  Geometrical  Theory  of 
Diffraction  (GTD) .  We  made  some  assumptions  like  the  scattering  body 
is  a  perfectly  conducting  body  surrounded  by  isotropic  and  homogeneous 
rtiedia,  the  radiations  are  monochromatic  and  have  high  frequency,  the 
dimensions  of  the  object  are  large  compared  to  the  wavelength,  and 
the  situation  is  static.  After  conducting  an  extensive  literature 
survey,  we  came  to  a  conclusion  that  GO  combined  with  GTD  is  the  most 
practical  way  of  calculating  scattering. 

We  applied  GO  and  GTD  to  calculate  backscatter  from  a  sphere 
and  from  a  spherical  segment  which  is  a  simple  representation  of  a 
car.  We  have  given  all  the  equations  in  detail  for  these  cases  and 
also  we  have  given  "Hierarchical  plus  Input- Process-Output"  charts 
for  the  ease  of  overview  and  coding  purposes.  These  examples  show 
that  GO  and  GTD  have  conceptual  simplicity  and  they  are  rather  easy 
to  apply. 

So  far  we  restricted  ourselves  to  the  scattering  from  perfectly 
conducting  surface  in  isotropic  and  homogeneous  media.  One  can  extend 
it  to  treat  partially  conducting  bodies  or  bodies  coated  with  dielect- 
tric  material  of  the  surrounding  being  inhomogeneous  and  anisotropic, 
or  any  combination  of  these.  Also,  there  is  a  lot  of  work  to  be  done 
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in  the  case  of  the  vertex  scattering.  It  is  one  of  the  challenging 
problems  to  find  the  diffraction  coefficient  for  a  vertex  or  tip 
diffraction.  One  can  study  the  effects  of  lower  frequency  on  the 
scattering.  Also  one  can  improve  upon  the  other  methods  of  the 
calculating  scattering  especially  the  numerical  technique  which  is 
more  promising. 


V. 
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APPENDIX  I 


HIERARCHY  PLUS  INPUT-PROCESS-OUTPUT 
(HIPO)  CHARTS 

HIPO  charts  facilitate  an  overview  of  the  system  and  determine 
the  program  functions  for  coding  purposes. 

The  HIPO  charts  that  follow  contain  a  visual  table  of 
contents  and  overview  diagrams  of  the  problem  that  has  been  dis¬ 


cussed. 
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